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Abstract 



The supersymmetric contributions to the Flavor Changing Neutral Current processes 
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may be suppressed by decoupling the scalars of the first and second generations. It is 
known, however, that the heavy scalars drive the stop mass squareds negative through the 
two-loop Renormalization Group evolution. This tension is studied in detail. Two new 
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^"'j items are included in this analysis: the effect of the top quark Yukawa coupling and the 

QCD corrections to the supersymmetric contributions to Am^- Even with Cabibbo-like 
degeneracy between the squarks of the first two generations, these squarks must be heavier 
than ~ 40 TeV to suppress Amy. This implies, in the case of a high scale of supersymmetry 
breaking, that the boundary value of the stop mass has to be greater than ~ 7 TeV to keep 
the stop mass squared positive at the weak scale. Low-energy supersymmetry breaking at 
a scale that is of the same order as the mass of the heavy scalars is also considered. In this 
case the finite parts of the two-loop diagrams are computed to estimate the contribution of 
the heavy scalar masses to the stop mass squared. It is found that for Cabibbo-like mixing 
between the squarks, the stop mass at the boundary needs to be larger than ~ 2 TeV. 
Thus, for both cases, the large boundary value of the stop masses leads to an unnatural 
amount of fine tuning to obtain the correct Z mass. 
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1 Introduction 



The origin of electroweak symmetry breaking (EWSB) and the subsequent gauge hierarchy prob- 
lem |l| are two large mysteries of the Standard Model (SM). Supersymmetry (SUSY)0 provides 
a promising solution to these problems, by both stabilising the weak scale against radiative 
corrections [0], and by naturally breaking the electroweak symmetry through the quantum cor- 
rections of the superpartner of the top quark to the Higgs boson mass ||. It is known, however, 
that generic weak scale values for the masses of the first two generation scalars give rates for 
many flavour violating processes that are in disagreement with the experimental observation. 
The measured value of Amx and detection limits for \i — > e-y, and fi — > 3e, for example, require 
that the first two generation scalars be degenerate to within a few tenths of a percent if their 
masses are at the weak scale |5], ||. Constraints from CP violation are generally even more severe. 
Understanding the origin of this degeneracy is the supersymmetric flavour problem. Attempts to 
resolve this puzzle without introducing any fine tuning include: using approximate non-abelian 
or abelian symmetries^]; communicating supersymmetry breaking to the visible sector by gauge 
interactions that do not distinguish between flavours |S]]; squark-quark mass matrix alignment 
H; and raising the soft masses of the first two generation scalars to the tens of TeV range 



lpj, y, y, y, y, y, y,|rj. 

The phenomenological viability and naturalness of this last scenario is the subject of this 
paper. To suppress flavour changing processes, the heavy scalars must have masses between a 
few TeV and a hundred TeV. The actual value depends on the degree of degeneracy and mixing 



between the masses of the first two generation scalars. As discussed in Reference [|18| , the masses 
of the heavy scalars cannot be made arbitrarily large without breaking colour and charge. This 
is because the heavy scalar masses contribute to the two-loop Renormalisation Group Equation 
(RGE) for the soft masses of the light scalars, such that the stop soft mass squared become more 
negative in RG scaling to smaller energy scales. This negative contribution is large if the scale 
at which supersymmetry breaking is communicated to the visible sector is close to the Grand 
Unification scale [ 18 1. With the first two generation soft scalar masses ~ 10 TeV, the initial value 
of the soft masses for the light stops must be ~ (few TeV) 2 to cancel this negative contribution 
|18| to obtain the correct vaccum. This requires, however, an unnatural amount of fine tuning 



to correctly break the electroweak symmetry [O, EU 



In this paper we analyse these issues and include two new items not previously discussed 
within this context: the effect of the large top quark Yukawa coupling, A t , in the RG evolution, 
that drives the stop soft mass squared more negative; and QCD radiative corrections in the 



Arrix constraint ||21|| . This modifies the bound on the heavy scalar masses which is consistent 
with the measured value of ArriK- This, in turn, affects the minimum value of the initial scalar 
masses that is required to keep the scalar masses positive at the weak scale. 

We note that the severe constraint obtained for the initial stop masses assumes that su- 
persymmetry breaking occurs at a high scale. This leaves open the possibility that requiring 
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positivity of the scalar masses is not a strong constraint if the scale of supersymmetry breaking 
is not much larger than the mass scale of the heavy scalars. In this paper we investigate this 
possibility by computing the finite parts of the same two-loop diagrams responsible for the neg- 
ative contribution to the light scalar RG equation, and use these results as an estimate of the 
two-loop contribution in an actual model of low-energy supersymmetry breaking. We find that 
in certain classes of models, requiring positivity of the soft masses may place strong necessary 
conditions that such models must satisfy in order to be phenomenologically viable. 

Our paper is organized as follows. In Section 2 an overview of the ingredients of our analysis 
is presented. Some philosophy and notation is discussed. Section 2.1 discusses the constraints 
on the masses and mixings of the first two generation scalars obtained from Amx after including 
QCD corrections. It is found, in particular, that Cabibbo-like mixing among both the first two 
generation left-handed squarks and right-handed squarks requires them to be heavier than 40 
TeV. Section 2.2 discusses the logic of our RG analysis, and some formulas are presented. This 
analysis is independent of the ArriK analysis. Sections 3 and 4 apply this machinery to the cases 
of low-energy and high-energy supersymmetry breaking, respectively. Section 3 deals with the 
case in which the scale at which SUSY breaking is communicated to the SM sparticles is close 
to the mass of the heavy scalars. We use the finite parts of the two-loop diagrams to estimate 
the negative contribution of the heavy scalars. We find that Cabibbo-like mixing among the 
left-handed and right-handed first two generation squarks implies that the boundary value of 
the stop masses has to greater than ~ 2 TeV to keep the stop mass squareds positive at the weak 
scale. This results in a fine tuning of naively 1% in electroweak symmetry breaking \I9\. We also 
discuss the cases where there is 0(1) mixing among only the right or left squarks of the first two 
generations, and find that requiring positivity of the slepton mass squareds implies a constraint 
on the stop masses of ~ 1 TeV if gauge-mediated boundary conditions are used to relate the 
two masses. This is comparable to the direct constraint on the initial stop masses. In Section 4, 
we consider the case where the SUSY breaking masses for the SM sparticles are generated at a 
high scale (~ 10 16 GeV). In this case, the negative contribution of the heavy scalars is enhanced 
by a large logarithm. We consider various boundary conditions for the stop and Higgs masses 
and find that with 0(0.22) degeneracy between the first two generation squarks, the boundary 
value of the stop mass needs to be larger than ~ 7 TeV. This gives a fine tuning of naively 



0.02%|19]. For 0(1) mixing between the left (right) squarks only, the minimum initial value of 
the stop is ~ 4(2) TeV. In Section 5 the scale of supersymmetry breaking is varied between 50 
TeV and 2 x 10 16 GeV. Uppers bounds on the amount of degneracy required between the first 
two generation scalars, that is consistent with positivity of the light scalar masses, naturalness 
in electroweak symmetry breaking, and the measured value of Am^-, are obtained. These results 
are summarized in Figures [12| and [H| We conclude in Section 6. In the Appendix, we discuss 
the computation of the two-loop diagrams which give the negative contribution of the heavy 
scalars to the light scalar mass squareds. 
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2 Overview. 



The chiral particle content of the Minimal Supersymmetric Standard Model (MSSM) contains 
3 generations of 5+10 representations of SU(5). The supersymmetry must be softly broken to 
not be excluded by experiment. Thus the theory must also be supplemented by some 'bare' soft 
supersymmetry breaking parameters, as well as a physical cutoff, M S usy- The 'bare' soft su- 
persymmetry breaking parameters are then the coefficients appearing in the Lagrangian, defined 
with a cutoff M S usy- It will be assumed for simplicity that the bare soft masses, fh^ , the bare 
gaugino masses Ma,o, and a bare trilinear term for the stops, AfA^o, are all generated close to 
this scale. The MSSM is then a good effective theory at energies below the scale Msusy, but 
above the mass of the heaviest superpartner. 

The physical observables at low-energies will depend on these parameters. If an unnatural 
degree of cancellation is required between the bare parameters of the theory to produce a mea- 
sured observable, the theory may be considered to be fine tuned. Of course, it is possible that a 
more fundamental theory may resolve in a natural manner the apparent fine tuning. The gauge- 
hierarchy problem is a well-known example of this. The Higgs boson mass of the SM is fine tuned 
if the SM is valid at energies above a few TeV. This fine tuning is removed if at energies close 
to the weak scale the SM is replaced by a more fundamental theory that is supersymmetric 0]. 

One quantification of the fine tuning of an observable O with respect to a bare parameter Aq 



is given by Barbieri-Giudice [19] to be 

A(0; A ) = (5O/O)/(5X /X ) = ^O. (1) 

It is argued that this only measures the sensitivity of O to Ao, and care should be taken when 
interpreting whether a large value of A necessarily implies that O is fine tuned [f2(J . It is not the 
intent of this paper to quantify fine tuning; rather, an estimate of the fine tuning is sufficient 
and Equation |l] will be used. In this paper the value of O is considered extremely unnatural if 
A(C; A ) > 100. 

The theoretical prediction for Am^ (within the MSSM) and its measured value are an exam- 
ple of such a fine tuning: Why should the masses of the first two generation scalars be degenerate 
to within 1 GeV, when their masses are 0(500 GeV)? Phrased differently, the first two gener- 
ation scalars must be extremely degenerate for the MSSM to not be excluded by the measured 
value of Am>K- An important direction in supersymmetry model building is aimed at attempting 
to explain the origin of this degeneracy. 

One proposed solution to avoid this fine tuning is to decouple the first two generation scalars 
since their masses are the most stringently constrained by the flavour violating processes JTTj 



IT, 12, [T^, |T|, [16|, [T7|]. In this scenario, some of the first two generation scalars have masses 



Ms 3> vtlz- To introduce some notation, n 5 (n w ) will denote the number of 5 (10) scalars of 
the MSSM particle content that are very heavy 0. Thus at energy scales E <C Ms the particle 
4 It is assumed that the heavy scalars form complete SU(5) multiplets to avoid a large Fayet-Illiopoulus D— 
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content is that of the MSSM, minus the n§ 5 and n w 10 scalars. In the literature this is often 
referred to as 'The More Minimal Supersymmetric Standard Model' [0. 

There are, however, other possible and equally valid sources of fine tunings. The measured 
value of the Z mass is such an example ||19|| . The minimum of the renormalized Higgs potential 



determines the value of the Z mass which is already known from experiment. The vacuum 
expectation value (vev) of the Higgs field is, in turn, a function of the bare parameters of the 
theory. The relation used here, valid at the tree-level, is 

2 mz ~ _/i + tan 2 /3-l • {2) 

It is clear from this Equation that requiring correct electroweak symmetry breaking relates the 
value of the soft Higgs masses at the weak scale, m 2 Hd (fj,G) and m^^fic), to the supersymmetric 
Higgs mass /i. A numerical computation determines the dependence of tt^hS^g) anc ^ m tf d (/ / G ! ) 
on the bare parameters Mao, m ^ n and Ms- In the MSSM, the cancellation required between 
the bare parameters of the theory for it not to be excluded by the Z mass increases as the 
scale of supersymmetry breaking is increased. Typically, the bare mass of the gluino, stops, 
and the first two generation squarks must be less than a few TeV and ten TeV, respectively, for 
successful electroweak symmetry breaking not to be fine tuned at more than the one per cent 
level H, fZg, 0]. 



These two potential fine tuning problems- the supersymmetric flavour problem and that of 
electroweak symmetry breaking- are not completely independent, for they both relate to the 



size of supersymmetry breaking |T7], Thus the consistency of any theoretical framework 

that attempts to resolve one fine tuning issue can be tested by requiring that it not reintroduce 
any comparable fine tunings in other sectors of the theory. This is the situation for the case 
under consideration here. Raising the masses of the first two generation scalars can resolve the 
supersymmetric flavour problem. As discussed in |L7|], this results in a fine tuning of m| through 



the two-loop dependence of rn 2 Hu (fiG) on M$. There is, however, another source of fine tuning 
of mz due to the heavy scalars: these large masses require that the bare masses of the stops, in 
particular, be typically larger than a few TeV to keep the soft masses squared positive at the 
weak scale |TB[. This large value for the bare stop mass prefers a large value for the vev of the 
Higgs field, thus introducing a fine tuning in the electroweak sector. Further, this fine tuning is 
typically not less than the original fine tuning in the flavour sector. This is the central issue of 
this paper. 

2.1 ArriK Constraints 

At the one-loop level the exchange of gluinos and squarks generates a AS = 2 operator. In the 
limit M 3 << Ms that we are interested in, the AS" = 2 effective Lagrangian at the scale Ms 



term at the one-loop level [p. 7, 13 
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obtained by integrating out the squarks is 

C eff = Y[^f ( Cl01 + + + C5O5 + h.c.) . (3) 

Terms that are 0(M^/Mg) are subdominant and neglected. We expand the exact result in 
powers of 5 L l,rr = s L . R c L , R r] L ^ R (fa{ - mf) LLfiR /fk 2 AVL R , where fh\ v is the average mass of the 
scalars, and where T)l,r is the phase and sl,r is the 1 — 2 element of the Wl,r matrix that appears 
at the gluino-squark-quark vertexf 5 ]. This approximation underestimates the magnitude of the 



exact result, so our analysis is conservative p8|| . The coefficients Cj to leading order in 6ll, $rr, 
are 

Ci = -225 d LL 
C A = 24S d LL S d RR 

C 5 = -A05 d LL 5 d RR . (4) 
The coefficient C\ is obtained from C\ with the replacement 5j jL — > 5 RR . The operators Oi are 

04 = d R SL, a d?L s R,b 

5 = d R s L , b dls R , a (5) 

and Oi is obtained from 0\ with the replacement L —>■ R. The Wilson coefficients, C\ — C5, are 
RG scaled from the scale of the squarks, Ms , to 900 MeV using the anomalous dimensions of the 



operators, 0\ — 0§. The anomalous dimension of 0\ is well known |22| and is fj,dC\/dfi = a s Ci/ir. 
We have computed the other anomalous dimensions and our result agrees with that of [^T|] (see 
this reference for a more general analysis of QCD corrections to the SUSY contributions to 
K — K mixing). These authors , however, choose to RG scale to Hhad, defined by a s (fihad)=^- 
The validity of the pertubation expansion is questionable at this scale; we choose instead to RG 
scale to 900 MeV, where a s (900 MeV) ~ .4. The result is 

C^had) = «lA(M 5 ) 

CtiUkad) = K A C 4 {M S ) + -(/C 4 - K 5 )C 5 (M S ) 

C^ had ) = k 5 C 5 (M s ) ' (6) 



where 



K4 = h 

.1/2 



/ x \ 6/27 / / x\6/25 / / x\6/23 / / n \ 6/21 / ,,, s \ 6/(9+(n 5 +3ni )/2) 

a s {m c ) \ / a s {m b ) \ / a s {m t ) \ / (Wg) \ / a s {M s ^ 
a s (900 MeV)/ \a s (m c ) J \a s (m b ) J \a s (m t ) J \a s (/j, c 

-4 



K 5 = k/ (7) 

5 In this paper only 1-2 generation mixing is considered. Direct L — R mass mixing is also neglected. 
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The effective Lagrangian at the hadronic scale is then 

Ceff = (-22(^ X )\ 1 C» 1 - 22(5 d RR )\ 1 d 1 + a^Of (4* + 5« 5 )C 4 - 40k 5 O 5 ) + h.c.) . 

(8) 

The SUSY contribution to the K — K mass difference is 

(Am K ) SUSY = 2Re < K\C eff \K > . (9) 
The relevant matrix elements (with bag factors set to 1) are 



<K\0 1 \K> = -m K f 



<K\O a \K> = + ± — m*/. 




p2 
'A' 



<K|0 5 |K> = + "7 )m K f K (10) 

\m s + m d J J 



in the vacuum insertion approximation. We use |23[ m^- = 497 MeV, = 160 MeV, m s = 150 
MeV, (Am K ) exp = 3.5x 10~ 12 MeV, anda s (M z ) = 0.118. This gives a s {m b ) = .21, a s (m c ) = .29 
and a s (900 MeV) = .38 using the one-loop RG evolution. Once values for (715, nio, Sf L , S RR ) are 
specified, a minimum value for Ms is gotten by requiring that (Ai71k)susy = (Am^exp- In the 
case that both 5 RR 7^ and 8ll 7^ 0, we assume that both the left-handed and right-handed 
squarks are heavy, so that (n 5 ,n 10 ) = (2,2). In this case we require that only the dominant 
contribution to Arrix, which is ~ 8f L 5 RR , equals the measured value of Attik- If 8rr 7^ and 
5ll = 0, we assume that only the right-handed squarks are heavy, and thus (n 5 ,ni ) = (2,0). 
Similarly, if 5ll 7^ and S RR = then (n 5 ,nio) — (0,2). Limits are given in Tables 1 and 2 for 



some choices of these parameters. These results agree with Reference [21 for the same choice 
of input parameters. For comparison, the limits gotten by neglecting the QCD corrections 
are also presented in Tables 1 and 2. We consider 5j jL (S RR ) = (i) 1, (ii) .22, (Hi) .1, and 
(iv) 0.04. These correspond to no mixing and no degeneracy; (ii) Cabibbo-like mixing; 
(Hi) Cabibbo-like mixing and ~ .5 degeneracy; and (iv) Cabibbo-like mixing and Cabibbo-like 
degeneracy. We expect only cases (i), (ii) and (Hi) to be relevant if the supersymmetric flavour 
problem is resolved by decoupling the first two generation scalars. From Table we note that 
for (n 5 , n w ) = (2, 0), Ms must be larger than ~ 30 TeV if it is assumed there is no small mixing 
or degeneracy (S RR = 1) between the first two generation scalars. 

The limits gotten from the mesaured rate of CP violation are now briefly discussed. Recall 
that the CP violating parameter e is approximately 



|Im < K\£ eff \K > 



(11) 



V2~Am K 

and its measured value is |e| ~ 1 7700 1 =2.3xl0 -3 [Q. In this case, the small value of e implies 
either that the phases appearing in the soft scalar mass matrix are extremely tiny, or that the 
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(n 5 ,ni ) = (2,2) 


(n 5 ,n w ) = (2,2) 




QCD incl. 


no QCD 


1 


182 TeV 


66 TeV 


.22 


40 TeV 


15 TeV 


.1 


18 TeV 


7.3 TeV 


0.04 


7.3 TeV 


3.1 TeV 



Table 1: Minimum values for heavy scalar masses Ms obtained from the measured value of Amx 
assuming M|/M| <C 1. The limits labeled 'QCD incl.' include QCD corrections as discussed in 
the text. Those labeled as 'no QCD' do not. 



MKr) (s d LL = o) 


(n 5 ,n w ) = (2,0) 


(n 5 ,n 10 ) = (2,0) 




QCD incl. 


no QCD 


i 


30 TeV 


38 TeV 


.22 


7.2 TeV 


8.9 TeV 


0.1 


3.4 TeV 


4.1 TeV 


0.04 


1.4 TeV 


1.7 TeV 



Table 2: Minimum values for heavy scalar masses M s obtained from the measured value of ArriK 
assuming Mf/Mj <C 1. The limits labeled as 'QCD incl.' include QCD corrections as discussed 
in the text. Those labeled as 'no QCD' do not. The limits for (n 5 ,n 10 ) = (0,2) obtained by 
5tt Srr are similar and not shown. 
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masses of the heavy scalars are larger than the limits given in Tables [1] and |2|. In the case 
where the phases are 0(1), Im < K\C e ff\K >~ Re < K\C e ff\K > and thus the stronger 
constraint on Ms is obtained from e and not Amjf, for the same choice of input parameters. 
In particular, the constraint from CP violation increases the minimum allowed value of Ms by 
a factor of \/\j2\/2t ~12.5. This significantly increases the minimum value of the initial light 
scalar masses that is allowed by the positivity requirement. 



2.2 RGE analysis 



The values of the soft masses at the weak scale are determined by the RG evolution. In the DR 
scheme p5|, [26|, ^7H , the RG equations []for the light scalar masses are, including the gaugino, A- 



term and \ t contributions at the one-loop level and the heavy scalar contribution at the two-loop 



level 28 



|m 2 (t = ln/i) = - ^ «a(^;m]W + ^ ^ Q«1W(«5^ + 3»iom? ) 

8 3 /4 3 1 \ 

+ T ^-(m% u (t) + mSg(t) + mj,(f) + A{tf) 

+ l^Y iai (t)TrYm 2 (t), (12) 
5 Air 

with r\ = (3, 2, 1) for = H u , t c , t, respectively, and zero otherwise. For simplicity it is assumed 
that Maa/ocaa are all equal at Msusy- The initial value of the gluino mass, M 3 , is then chosen 
to be the independent parameter. To avoid a large Fayet-Illiopoulus D-term at the one-loop 
level, we assume that the heavy scalars form complete SU(5) representations |17, [R|. We use 



SU(5) normalisation for the U(l) coupling constant and Q = T 3 + Y. Finally, C\ is the quadratic 
Caismir for the gauge group Ga that is 4/3 and 3/4 for the fundamental representations of 577(3) 
and 577(2), and 3/5Y? for the U(l) group. The cases (n 5 ,n 10 )= (I) (2,2), (II) (2,0), (III) (0,2) 
are considered. The results for the case (3, 0) is obtained, to a good approximation, from Case 
(II) by a simple scaling, and it is not discussed any further. 

Inspection of Equation ([E|) reveals that in RG scaling from a high scale to a smaller scale the 
two-loop gauge contribution to the soft masses is negative, and that of the gauginos is positive. 
The presence of the large \ t Yukawa coupling in the RGE drives the value of the stop soft mass 
squared even more negative. This effect increases the bound on the initial value for the stop 
soft masses and is included in our analysis. In our analysis the top quark mass in MS scheme 
is fixed at 167 GeV. 



6 An earlier version of the analysis presented in this manuscript did not include the one-loop hypercharge 
I?— term in the RG equations. This has been corrected. This has not changed our conclusions though, since the 
numerical effect of this error is small. The authors thank B. Nelson for drawing our attention to this omission. 
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In the MSSM there is an extra parameter, tan/?, which is the ratio of the vacuum expeca- 
tions values of the Higgs fields that couple to the up-type and down-type quarks respectively. 
Electroweak symmetry breaking then determines the top quark mass to be m t = \ t / \f2v sin (3 
with v ~ 247 GeV. In our analysis we consider the regime of small to moderate tan (3, so that all 
Yukawa couplings other than At are neglected in the RG evolution. In this approximation the 
numerical results for factor t c are independent of tan f3. In the numerical analysis of Sections 
|3] and | tan/5=2.2 is considered. In Section |]tan/3 = 10 is also considered. 

In the case of low-energy supersymmetry breaking, the scale Msusy is not much larger than 
the mass scale of the heavy scalars. Then the logarithm ~ln (Msusy /Ms) that appears in the 
solution to the previous RG equations is only 0(1). In this case the finite parts of the two-loop 
diagrams may not be negligible and should be included in our analysis. We use these finite parts 
to estimate the size of the two-loop heavy scalar contribution in an actual model. 

The full-two loop expression for the soft scalar mass at a renormalisation scale hr is tn 2 ull (jj, R ) = 
m|^/(/iR) + m 2 j inite (nR}, where is the solution to the RG equation in DR scheme, and 

m 'finite(^R) * s the finhe part of the one- loop and two- loop diagrams, also computed in DR 
scheme. The finite parts of the two-loop diagrams that contain internal heavy scalars are com- 
puted in the Appendix and the details are given therein. The answer for these two-loop finite 
parts is (assuming all heavy scalars are degenerate with common mass M|) 

ml fintt M = -i(ln(47r) - 7 + y - 2 - In (^|)) x £ (^^Y(n 8 + 3n 10 )^M| 

3 1 ( 2 A/f^ 

~5 16^ aiW(n5 ~ nio)Yi ( 6 ~ 3 71 " 2 + 2 ( ln ( 47r ) " i) ~ 41n( ^f ) 

x fgOs^fl) - -a 2 {fi R ) - M 2 S (13) 

where the gaugino and fermion masses are neglected. Since we use the DR' scheme to compute 
the finite parts of the soft scalar masses, the limits we obtain on the initial masses are only 
valid, strictly speaking, in this scheme. This is especially relevant for the case of low scale SUSY 
breaking. So while these finite parts should be viewed as semi-quantitative, they should suffice 
for a discussion of the fine tuning that results from the limit on the bare stop mass. For the 
case of high scale SUSY breaking, the RG logarithm is large and so the finite parts are not that 
important. 

Our numerical analysis for either low-energy or high-energy supersymmetry breaking is de- 
scribed as follows. 

The RG equations are evolved from the scale Msusy to the scale at which the heavy scalars 
are decoupled. This scale is denoted by us and should be 0(M S ). The RG scaling of the heavy 
scalars is neglected. At this scale the finite parts of the two-loop diagrams are added to m% (/is). 

Ji 

We note that since the two-loop information included in our RG analysis is the leading O(Mj) 
effect, it is sufficient to only use tree-level matching at the scale us- Since the heavy scalars are 
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not included in the effective theory below Ms and do not contribute to the gauge coupling beta 
functions, the numerical results contain an implicit dependence on the number of heavy scalars. 
This results in a smaller value for 0:3 (/is) compared to its value if instead all the scalars have a 
~ ITeV mass. This tends to weaken the constraint, and so it is included in our analysis f\ The 
soft masses are then evolved using the one-loop RGE to the mass scale at which the gluinos are 
decoupled. This scale is fixed to be /xg=1 TeV. 

A constraint on the initial value of the soft masses is obtained by requiring that at the weak 
scale the physical scalar masses are positive. The experimental limit is ~ 70 GeV for charged 
or coloured scalars [^J. The physical mass of a scalar is equal to the sum of the soft scalar 
mass, the electro- weak D— term, the supersymmetric contribution, and some finite one-loop 
and two-loop contributions. As mentioned in the previous paragraph, in the effective theory 
below Ms the finite two-loop part from the heavy scalars is included in value of the soft scalar 
mass of the light sparticles at the boundary, defined at \ir = us ~ M s . The finite one-loop 
contributions are proportional to the gaugino and other light scalar masses, and are smaller 
than the corresponding logarithm that is summed in m 2 (/i/j). So we neglect these finite one- 
loop parts. Further, the electroweak D— terms are less than 70 GeV. For the scalars other than 
the stops, the supersymmetric contribution is negligible. In what follows then, we will require 
that fhj(fj,G) > for scalars other than the stops. The discussion with the stops is complicated 
by both the large supersymmetric contribution, m 2 , to the physical mass and by the L — R 
mixing between the gauge eigenstates. This mixing results in a state with mass squared less 
than min(m~ + m 2 , m~ c + m 2 ), so it is a conservative assumption to require that for both gauge 
eigenstates the value of m~ +m? is larger than the experimental limit. This implies that m 2 ^(70 
GeV) 2 -(175 GeV) 2 = -(160 GeV) 2 . In what follows we require instead that m\ > 0. This 
results in an error that is (160GeV) 2 /2mj. o m 26 GeV if the constraint obtained by neglecting 
rrit is ~ 1 TeV. For the parameter range of interest it will be shown that the limit on the initial 
squark masses is ~ 1 TeV, so this approximation is consistent. 

We then combine the above two analyses as follows. The Am^- constraints of Section 12" 



determine a minimum value for Ms once some theoretical preference for the 5's is given. Either a 
natural value for the 5's is predicted by some model, or the 5's are arbitrary and chosen solely by 
naturalness considerations. Namely, in the latter case the fine tuning to suppress Amx is roughly 
2/5. Further, a model may also predict the ratio M3/M5. Otherwise, Equations [l] and ^| may be 
used as a rough guide to determine an upper value for M 3 , based upon naturalness considerations 
of the Z mass. Without such a limitation, the positivity requirements are completely irrelevant 
if the bare gluino mass is suffuciently large; but then the Z mass is fine tuned. Using these 
values of M3 and Ms, the RGE analysis gives a minimum value for the initial stop masses which 
is consistent with Amx and positivity of the soft masses. This translates into some fine tuning 



7 This is the origin of a small numerical discrepancy of ~ 10% between our results and the analysis of |18[ in 
the approximation A t = 0. 
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of the Z mass, which is then roughly quantified by Equations [T] and |2|. 

Finally, we remark that our analysis may also be extended to include models that contain a 
Fayet-Illiopoulos hypercharge D— term, Cd-, at the tree-level. The effect of the D— term is to shift 
the soft scalar masses, m 2 — > m 2 = m 2 + Yi( D . In this case, the positivity analysis applies to 
ra 2 , rather than m 2 . 

3 Low Energy Supersymmetry Breaking 

In this Section we investigate the positivity requirement within a framework that satisifes both of 
the following: (i) supersymmetry breaking is communicated to the visible sector at low energies; 
and (ii) multi-TeV scale soft masses, Ms, are generated for some of the first two generation 
scalars. This differs from the usual low-energy supersymmetry breaking scenario in that we 
assume M| 3> m~ Q . In the absence of a specific model, however, it is difficult to obtain from 
the positivity criterion robust constraints on the scalar spectra for the following reasons. At 
the scale Msusy it is expected that, in addition to the heavy scalars of the MSSM, there are 
particles that may have SM quantum numbers and supersymmetry breaking mass parameters. 
All these extra states contribute to the soft scalar masses of the light particles. The sign of this 
contribution depends on, among other things, whether the soft mass squared for these additional 
particles is positive or negative-clearly very model-dependent. The total two-loop contribution 
to the light scalar masses is thus a sum of a model-dependent part and a model independent part. 
By considering only the model-independent contribution we have only isolated one particular 
contribution to the total value of the soft scalar masses near the supersymmetry breaking scale. 
We will, however, use these results to estimate the typical size of the finite parts in an actual 
model. That is, if in an actual model the sign of the finite parts is negative and its size is of the 
same magnitude as in Equation (|H|), the constraint in that model is identical to the constraint 
that we obtain. The constraint for other values for the finite parts is then obtained from our 
results by a simple scaling. 

Before discussing the numerical results, the size of the finite contributions are estimated in 
order to illustrate the problem. Substituting M s ~ 25 TeV, 013(25 TeV) ~ 0.07 and ai(25 
TeV) ~ 0.018 into Equation [13| gives 

/ M \ 2 

5m\ » -(410 GeV) 2 (n 5 + 3n 10 ) (— |y) (14) 

for squarks, and 

25tW) (15) 

for the right-handed selectron. The negative contribution is large if M5 ~ 25 TeV. For example, 
if n 5 = mo = 2 then 5m% « -(200 GeV) 2 and bm\ « -(1.2 TeV) 2 . If n 5 = 2, n w = 0, then 
5ml ~ -(170 GeV) 2 and 5m\ w -(580 GeV) 2 . 
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In this low-energy supersymmetry breaking scenario, it is expected that Msusy ~ Ms- In our 
numerical analysis we will set Msusy = A*s since the actual messenger scale is not known. The 
scale fis is chosen to be 50 TeV. At the scale /is =50 TeV the /is-independent parts of Equation 
(|13|) are added to the initial value of the soft scalar masses. The soft masses are then evolved 
using the RG equations (not including the two-loop contribution) to the scale hq— ITeV. 

First we discuss the constraints the positivity requirement imply for /; ^ or £r. In this 
case m| is renormalised by M| , Mj, m 2 j and the initial value of TrFm 2 = Dy.o- We find 

m|0 G ) = m| + (0.243C| + 0.0168Q + 0.00156lf )M 3 2 + c D x lO^L^o 

-(0.468C£ + .095C2 + .0173lf )i(n 5 + 3n 1Q ) x 10~ 3 M| 
-0.0174(n 5 - n w )Yi x 10" 3 Mj 

-(n 5 - mo) ((-0.00058 + 0.0016(n 5 + 3m ))Mf - .925M 3 2 ) x 10" 3 , (16) 

where the strongest dependence on (n 5 , nio) has been isolated. The coefficient appearing in front 
of -Dy,o is cd = —6. The numerical coefficients in Equation ( |16|) also depend on (n5,nio) and the 



numbers presented in Equation (jig) are for (ns,nio) = (2,0). This sensitivity is, however, only a 
few percent between the four cases under consideration here Requiring positivity of the soft 
scalar masses directly constrains m 2 - /M s and M| /M|. 

The value of _Dy,o depends on the spectrum at the supersymmetry breaking scale, and is 
therefore model-dependent. To obtain model-independent constraints from the positivity re- 
quirement, we therefore only constrain the combination m| = m| + CoYi-Dyo- Only this 
combination appears in the weak-scale value for the scalar mass of fi. The numerical effect is 
small, since with Dy.o ~ 0(m 2 f ), the coefficient of m| is shifted from 1 to ~ 1 — 6. x 10 -3 Yj. 

The positivity requirement m| for fi 7^ t or t c is given in Figure p] for different values of 

Ji 

n 5 and n w . That is, in Figure [l] the minimum value of m^ /M s required to keep the soft 
masses positive at the scale hg is plotted versus M 3t0 /Ms- We conclude from these figures that 
the positivity criterion is weakest for n§=2 and nio=0. This is expected since in this case the 
heavy particle content is the smallest. We note that even in this 'most minimal' scenario the 
negative contribution to the masses are rather large. In particular, we infer from Figure [I] that 
for (n 5 = 2, mo = 0) and M s ~ 25 TeV, 5m% ss -(190 GeV) 2 for M 3i0 as large as 1 TeV. In 
this case it is the two-loop contribution from the hypercharge D-term that is responsible for the 
large negative mass squared. In the case (71,5, nio) = (2, 2), we obtain from Figure [I] that for Ms ~ 
25 TeV, 5m% w -(210 GeV) 2 and 5m\ c w -(1.1 TeV) 2 for M 3i0 as large as 1 TeV. 

We now apply the positivity requirement to the stop sector. In this case it is not possible to 
directly constrain the boundary values of the stops for the following simple reason. There are 
only two positivity constraints, whereas the values of m 2 (/xc) and m 2 c (fiG) ar e functions of the 
three soft scalar masses fn 2 , m~ c Q and rn 2 Hu . To obtain a limit some theoretical assumptions 
8 This dependence is included in Figure |l|. 
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Figure 1: Limits for mi Q /M$ from the requirement that the mass squareds are positive at the 
weak scale, for low-energy supersymmetry breaking. The regions below the curves are excluded. 
For the case (2,0), the limits for the other squarks are very similar to that for Q and are therefore 
not shown. 

must be made to relate the three initial soft scalar masses. 

The numerical solutions to the RG equations for tan/3=2.2 and (n 5 ,nio) = (2,0) are 

mf(> G ) = -0.0303A 2 + 0.00997AM 3i0 + 0.322M 3 2 + c D x - x 10" 3 Dy i0 

-0.0399(m^ uj0 + m% fi ) + 0.960m 2 - 0.000645c L M| 

_2 

m%(u G ) = -0.0606A 2 + 0.0199A t M 3i0 + 0.296M 3 2 + c D x — x lO" 3 ^ 
0.920m 2 C;0 - 0.0797(m^ + m| ) - 0.000492c R M| 

m 2 Hu (u G ) = -0.0909 A 2 + 0.0299 A t M 3fi - 0.0289M 2 + c D x - x 10~ 3 LV,o 

+0.880m| u - 0.119(m| + mf co ) + 0.0000719c H M|. (17) 

The numerical coefficients other than that of Ms do not vary more than a few percent between 
the different values for (n^,,nio), and thus this dependence is not shown. For Ms, we find that 
(c L ,c R ,c H ) is (1,1,1), (3.62,3.84,4.59), (2.78,3.04,3.92), for (n 5 ,n 10 ) = (2,0), (2,2) and (0,2), 
respectively. Also, = —6. We find from Equations |TJ and |2| that to keep m 2 z fine tuned at less 
than 1% (A < 100) in each of the bare parameters, we must have: 460 GeV; M 3:0 ^2.3 TeV; 
mj ^1.7 TeV; m5^80 TeV and mi ^50 TeV for (ns,nio) = (2,2). Finally, for other values of 
these parameters the fine tuning increases as A = 100 x m 2 fm\, where m Q is the value of in that 
gives A = 100. 
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It is possible to show, using the fact that Yh u + Yq + Y u c = 0, that the solutions in Equations 
17] are unchanged if we replace m 2 with m 2 = m 2 + CD X 10~ 3 Y;ZV,o and set _Dy,o — 0. In what 



follows then, we will use the posivitity analysis to constrain rh 2 for the stops. We note though, 
that the difference between m 2 Q and m 2 is small, owing to the small coefficient appearing in 
front of -Dy,o- In the remainder of this Section the tilde on m 2 will be removed to simplify the 
notation. 

To constrain the initial values of the stop masses we will only consider gauge-mediated su- 
persymmetry breaking mass relations. From Equation [L7] we see that to naturally break elec- 
troweak symmetry a small hierarchy m~ > rn 2 Hu is required. This is naturally provided by 
gauge-mediated boundary conditions f\ The relations between the soft scalar masses when su- 
persymmetry breaking is communicated to the visible sector by gauge messengers are || 



3 v t a 2 A (M SUSY ) 

4^ A a 2 {M SUSY ) + a 2 {M SUSY )/b tc 



Substituting these relations into Equations (|T7D and assuming A tt o =0 determines m 2 (nc) and 
^^(/^g) as a function of M^o, Mj and m~ c Q . In Figure we have plotted the minimum value of 
m t - c o /M 3i o required to maintain both m|(/ic) > and m~ c (/iG) > 0. 

Another interesting constraint on these class of models is found if it is assumed that the initial 
masses of all the light fields are related at the supersymmetry breaking scale by some gauge- 
mediated supersymmetry breaking (GMSB) mass relations, as in Equation fll8|) . This ensures 
the degeneracy, as required by the flavour changing constraints, of any light scalars of the first 
two generations. This is required if, for example, one of n 5 or n w are zero. Then in our previous 
limits of mj. for fa ^ t or t c , constraints on the initial value of m f - c are obtained by relating mj. 
to m£c o using Equation (|18|). In this case the slepton masses provide the strongest constraint and 
they are also shown in Figure 0. This result may be understood from the following considerations. 
The two-loop hypercharge D-term contribution to the soft mass is ~ Yi(n 5 — n^aia^Mg and 
this has two interesting consequences. The first is that for ^ nio, the resulting 5m 2 is always 
negative for one of e c or L. Thus in this case there is always a constraint on m~ c once gauge- 
mediated boundary conditions are assumed. That this negative contribution is large is seen 
as follows. The combined tree-level mass and two-loop contribution to the selectron mass is 
approximately m% — kaia 3 Mg where k is a numerical factor. Substituting the gauge-mediated 
relation m 2 C Q ~ a\/a\m 2 c , the combined selectron mass is al/al(m 2 c Q — k{a 3 / ai)a\Mg) . Since 
the combined mass of the stop is ~ m 2 c Q — fc'a|Mj, the limit for m~ c obtained from the positivity 
requirement for m 2 c is comparable or larger than the constraint obtained from requiring that 
m~ c remains positive. For example, with n 5 = 2, n w = and Ms ~ 25 TeV, the right-handed 
slepton constraint requires that m t ~ c ~ 1.1 TeV. For n 1Q =2, n 5 =0 and M s ~ 25 TeV, L 3 is 

9 In fact, low-energy gauge-mediated supersymmetry breaking provides "too much" electroweak symmetry 
breaking p9|. 
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Figure 2: Limits for mic^/Ms from the requirement that the stop and slepton mass squared 
are positive at the weak scale. The regions below the curves are excluded. Low-energy gauge- 
mediated supersymmetry breaking mass relations between the light sparticles and tan/3 =2.2 
are assumed. 

driven negative and implies that mjc ~ 1 TeV. From Figure |2| we find that these results are 
comparable to the direct constraint on m t - c obtained by requiring that colour is not broken. 

The positivity analysis only constrains mi. o/Ms for a fixed value of M^^/Ms- To directly 
limit the initial scalar masses some additional information is needed. This is provided by the 
measured value of Amx- If some mixing and degeneracy between the first two generation scalars 
is assumed, parameterized by (8ll,5rr), a minimum value for Ms is obtained by requiring 
that the supersymmetric contribution to Am^ does not exceed the measured value. We use 
the results given in Section ^| to calculate this minimum value. This result together with the 
positivity analysis then determines a minimum value for m 4 - c for a given initial gluino mass 
M 3 . The RG analysis is repeated with /is = M s , rather than /i 5 =50 TeV. We only present the 
results found by assuming GMSB mass relations between the scalars. These results are shown 
in Figure |3|. The mass limits for other fa are easily obtained from the information provided 
in Figure [I] and Table |] and are not shown. From Figure ||] we find that for (n.5, nio) = (2, 2) 
and M 3:0 less than 2 TeV, mf c0 must be larger than 8 TeV for \/^ljJrr = 1, and larger than 
1.8 TeV for v^Zl^rr = -22. This results in A(m|,m~ ) of 2000 and 120, respectively In this 
case both the squark and selectron limits for comparable. The limits for other choices 

for \/5ll5rr are obtained from Figure |3] by a simple scaling, since to a good approximation 
Arrix ~ 5ll&rr/ M s . For the cases (ri5,n w ) = (2,0) and (0,2), the corresponding limits are 
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much weaker. In the case (n 5 ,n w ) = (2,0), for example, only for 5rr ~ 1 does the selectron 
mass limit require that m t - c ~1 TeV. The limits for a smaller value of 6 are not shown. 

We conclude with some comments about how these results change if CP violation is present 
in these theories with 0(1) phases. Recall from Section that for the same choice of input 
parameters, the limits on the initial stop masses increases by about a factor of 12. This may 
be interpreted in one of two ways. Firstly, this constrains those models that were relatively 
unconstrained by the Aitlk limit. We concentrate on those models with = 2 and n\o = 0, 
since this case is the most weakly constrained by the combined Amx and positivity analysis. 
The conclusions for other models will be qualitatively the same. We find from Figure |3] the limit 
mj c o >1 TeV 0is only true if 5rr ~ 0(1). Smaller values of Srr do not require large initial stop 
masses. From the CP violation constraint, however, smaller values for Srr are now constrained. 
For example, if 5rr ~0.1 and 0(1) phases are present, then mje >1 TeV is required. Secondly, 
the strong constraint from e could partially or completely compensate a weakened constraint 
from the positivity analysis. This could occur, for example, if in an actual model the negative 
two-loop contribution to the stop mass squared for the same initial input parameters is smaller 
than the estimate used here. For example, if the estimate of the two-loop contribution in an 
actual model decreases by a factor of ~ (12. 5) 2 and 0(1) phases are present, the limit in this 
case from e for the same 5 is identical to the values presented in Figure |3|. 



4 High Scale Supersymmetry Breaking 

In this section, we consider the case in which SUSY breaking is communicated to the MSSM 
fields at a high energy scale, that is taken to be Mqut = 2 x 10 16 GeV. In this case, the 
negative contribution of the heavy scalar soft masses to the soft mass squareds of the light 
scalars is enhanced by ~ \n(M GUT /50 TeV), since the heavy scalar soft masses contribute to the 
RGE from Mqvt to mass of the heavy scalars. It is clear that as the scale of SUSY breaking is 
lowered the negative contribution of the heavy scalar soft masses reduces. 

This scenario was investigated in Reference ]18[ , and we briefly discuss the difference between 



that analysis and the results presented here. In the analysis of Reference ||18|| , the authors made 
the conservative choice of neglecting X t in the RG evolution. The large value of At can change 
the analysis, and it is included here. We find that for some pattern of initial stop and up-type 
Higgs scalar masses, e.g. universal scalar masses, this effect increases the constraint on the stop 
masses by almost a factor of two. This results in an increase of a factor of 3-4 in the amount of 
fine tuning required to obtain the correct Z mass. Further, in combining the positivity analysis 



10 For GMSB relations only. The direct constraint on the stop masses is slightly weaker. 

11 This choice for the high scale is done to remain agnostic about any physics appearing between the Grand 
Unification scale and the Planck scale. This also results in a conservative assumption, since the negative two-loop 
contribution is smaller with Msusy = Mqut- 
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Figure 3: Limits for m t ~ c from the requirement that the stop and slepton mass squared are 
positive at the weak scale while suppressing ArriK, for different values of (715, nio), and (5ll, 8rr)- 
The regions below the curves are excluded. Low-energy gauge-mediated supersymmetry breaking 
mass relations between the light scalars and tan/? =2.2 are assumed. 



17 



with the constraints from the Amx analysis, the QCD corrections to the Flavour Changing 
Neutral Current (FCNC) operators has been included, as discussed in Section |2|. In the case 
(n5,n 10 ) = (2,2), this effect alone increases the positivity limit by a factor of ~ 2 — 3. The 
combination of these two elements imply that the positivity constraints can be quite severe. 

We proceed as follows. First, we solve the RGEs from Mqut to fis where the heavy scalars 
are decoupled. At this scale, we add the finite parts of the two-loop diagrams. Next, we RG 
scale (without the heavy scalar terms in the RGEs) from fis to p,Q using these new boundary 
conditions. Except where stated otherwise, the scales fig and p,Q are fixed to be 50 TeV and 1 
TeV, respectively. 

For fi 7^ t, t c we find, 

m | (^ G ) = m| Q + (2.84C* + 0.639Q + 0.159>f )M 3 2 + c D YiD Yfi 

-(4.38C£ + 1.92C* + 0.622lf )-(n 5 + 3n 10 ) x 10~ 3 M| 
-0.829(t7 5 -n w )Yi x 10~ 3 M| 

+{n 5 - mo) (l7.2M 3 2 + (.226 - 0.011(n 5 + 3n 10 ))M|) Y x 10~ 3 . (19) 



These results agree with Reference JL8| for the same choice of input parameters. The term pro- 
portional to Dy,o, and the terms in the last line result from integrating the one-loop hypercharge 
D— term. In this case cp = —0.051. As in the previous Section, the numerical coefficients in 
Equation ([19]) depend on (n 5 , n w ) through the gauge coupling evolution, and the numbers in 
Equation(|T9"D are for (n 5 ,n w ) = (2,0) 0. Requiring the soft masses squared to be positive 
constrains rhf = mf + cd^-D^o- In Figure § we plot the values of rh^ io /Ms that determine 
m~ (hq) = as a function of M 3 /M s , for fa = L iy Q iy d\ and e\. We emphasize that the 

Ji 

results presented in Figure f| are independent of any further limits that FCNC or fine tuning 
considerations may imply, and are thus useful constraints on any model building attempts. 
For the stops, the numerical solutions to the RGEs for tan f3 = 2.2 are 

m 2 (/i G ) = -0.021A 2 + 0.068A t M 3 , + 3.52M 3 2 + c D ^D Y0 

-0.142(77^ + mf C)0 ) + 0.858m? - c L 0.00613M| 

m 2 c (/i G ) = -0.042A 2 + 0.137A t M 3 , + 2.33M 3 2 + c D ^D Y0 
-0.283(777^ + m| ) + 0.716m| - c R 0.00252M| 

mlr» = -0.063A 2 + 0.206AM3,o-l-72M 3 2 + c D ^ y>0 

-0.425(mf o + mf Cj0 ) + 0.574777| u o + c H 0.00193M| (20) 

where (c L ,c R ,c H ) = (1,1,1), (3.57,4.92,5.15), (2.7,4.16,4.27) for (r7 5 ,77 10 ) = (2,0), (2,2) and 
(0,2), respectively. Also, cd = —0.051. The mixed two-loop contribution to the RG evolution is 
12 The numerical results presented in Figure |4| include this dependence. 
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Figure 4: Limits for mj. for different values of (ns,nio) from the requirement that the mass 
squareds are positive at the weak scale, assuming a supersymmetry breaking scale of Mqut- The 
regions below the lines are excluded. 

oc (715 — nio) and is not negligible. Thus there is no simple relation between the c's for different 
values of and n\Q. From Equations [2] and |1] we find that to keep m 2 z fine tuned at less than 
1% (A < 100) in each of the bare parameters, we must have: /x^ 460 GeV; M 30 ^300 GeV; 
m ( - ^.87 TeV; m 5 j^l6 TeV; and m 10 ^10 TeV, for (n 5 ,n 10 ) = (2,2). The fine tuning of the 
Z mass with respect to the heavy scalars is discussed in |T7[. Finally, for other values of these 
parameters the fine tuning increases as A = 100 x rh 2 /rfiQ, where m is the value of m that gives 
A = 100. 

As in Section RL we rewrite Equations |2D| in terms of m 2 Q = m 2 Q + CDYiD Y fl- This is equivalent 
to setting D Y ,o = in Equations [2(], and relabeling m 2 — > m 2 . In what follows, we use the 
positivity analysis to constrain m 2 Q . Since cd is small and Dy,o ~ 0(m 2 ), the difference between 
rh 2 and m 2 is small. To simplify the notation, in the remainder of this Section we will also 
remove the tilde from fh 2 . 

As was also discussed in Section [|, some relations between m~ , m 2 c and rn 2 Hu are needed 
to obtain a constraint from Equation(pO|), using to|(//g) > and m| c (/xc) > 0. We discuss both 
model-dependent and model-independent constraints on the initial values of the stop masses. The 
outline of the rest of this Section is as follows. First, we assume universal boundary conditions. 
These results are presented in Figure Model-independent constraints are obtained by the 
following. We assume that m^ H = and choose A t)0 to maximize the value of the stop masses 
at the weak scale. These results are presented in Figure ^. We further argue that these constraints 
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represent minimum constraints as long as rn 2 Hu > 0. To obtain another set of model independent 
constraints, we use the electroweak symmetry breaking relation to eliminate m 2 Hu in favour of 
fx. Then we present the positivity limits for different values of jl/Mg, where jl 2 = fi 2 + \rn? z , 
and assume that m 2 ld = to minimize the value of fi 0. These limits are model-independent 
and are presented in Figure 0, for the case = uiq = 2. We then combine these analyses with 
the limits on Ms obtained from Am^. We conclude with some discussion about the anomalous 
D— term solutions to the flavour problem. 

We first consider universal boundary conditions for the stop and Higgs masses. That is, we 
assume that m~ = m~ c Q = m 2 Hu = m^. In Figure [5] we plot for tan j3 = 2.2 the minimum value 
of rho/Ms required to maintain m 2 (/ic) > and m~ a {fiG) > 0. This value of tan (3 corresponds to 
Xt(McuT) = -88, in the case that (n 5 , n w ) = (2, 0). For comparison, the results gotten assuming 
X t = may be found in Reference [|18[]. For n§ = Uiq = 2 we note from Figure [5] that if Ms = 20 
TeV and the gaugino masses are small, the limit on the stop mass is m~ tc > 6.2 TeV. This limit 
is weakened to 6 TeV if M 3j0 ^ 300 GeV is allowed. Even in this case, this large initial stop mass 
requires a fine tuning that in this case is A ~ (6 TeV) 2 /m| ~ 4200, i.e. a fine tuning of ^10~ 3 
is needed to obtain the correct Z mass. 

We now assume rn 2 Hu = and choose the initial value of A t $ to maximize the value of 
m~ . (/ic). The values of m~ Q and chosen such that m?(/ic) > and m 2 c (fj,G) > 0. We 

note that in this case the constraint is weaker because the Xt contribution to the RG evolution 
of the stop masses is less negative. These results are plotted in Figure ||. 

We discuss this case in some more detail and argue that the minimum value of obtained 
in this way will be valid for all m 2 Hu > and all A t0 . Eliminate the A t term by choosing 
A t fl = kMzfl such that the A t contributions to m~ . (a*g) is maximized. Other choices for A t £ 
require larger values for m 2 Q to maintain m 2 (fie) — 0. The value of k is determined by the 
following. A general expression for the value of the soft masses of the stops at the weak scale is 



™$(hg) = -aA 2 o + bA tfi M 3fi + cM 2 fi + ■■■, (21) 

m%(jiG) = -2a^ + 2bA tfi M 3fi + dMl fi + ■ ■ • , (22) 

with a, c and d positive. The maximum value of m~ (hg) is obtained by choosing A t ^ = hM^fl/2a. 
The value of the stops masses at this choice of A t are 

m?(^) = (c+^)M 3 2 + ---, (23) 
m|(/i G ) = (rf + 2^)M3 2 + ---. (24) 



An inspection of Equation 20 gives b = 0.068 and a = 0.021 for tan/3 = 2.2. In this case the 



'best' value for A t $ is Af ~ 1.6M 30 . It then follows that the quantity b 2 /4a = 0.055 is a small 

13 Strictly speaking, this last assumption is unnecessary. Only the combination jl 2 H = jl 2 — m 2 Hd / tan 2 (3 appears 
in our analysis. Thus for m 2 Hd 7^ our results are unchanged if the replacement jl — > jlu is made. 
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correction to the coefficient of the gaugino contribution in Equation 20. Thus the difference 
between the minimum initial stop masses for A ti o = and A t> o= Af is small. Next assume that 
m 2 Hu o = 0. Requiring that both m?(/ic) = and m~ c (//(j) = determines a minimum value for 



ni 



j and m 2 CQ . Now since the rn 2 Hu contribution to both the stop soft masses is negative (see 
Equation ^0|), the minimum values for m| found by the preceeding procedure are also minimum 
values if we now allow any m 2 Hu > 0. 

We conclude that for all A t $ and all rn 2 Hu > 0, the limits presented in Figure |6| represent 
lower limits on the initial stop masses if we require that the soft masses remain positive at the 
weak scale. Further, the limits in this case are quite strong. For example, from Figure |6| we 
find that if M s ~ 20 TeV and M 3i0 ~ 200 GeV (so that M 3fi /M s ~10" 2 ), then the initial stop 
masses must be greater than 3.5 TeV in the case that (n 5 ,n w ) = (2,2) The results are stronger 
in a more realistic scenario, i.e. m 2 Hu0 > 0. If, for example, m 2 Hu0 ~ m 2 c0 /9 the constraints 
are larger by only a few percent. In the case that , m 2 Hufi = m 2 c Q = rri~ Q , presented in Figure |], 
however, the constraint on the initial t c mass increases by almost a factor of two. 

To obtain constraints on the initial stop masses we have thus far had to assume some relation 
between m 2 Hu and m~ c ; e.g., rn 2 Hu = or rn 2 Hu Q = m~ c Q . Perhaps a better approach is to use 
the EWSB relation, Equation (^]), to eliminate rn 2 lu0 in favour of /i 2 . This has the advantage 
of being model- independent. It is also a useful reorganization of independent parameters since 
the amount of fine tuning required to obtain the correct Z mass increases as \i is increased. To 
obtain some limits we choose rn 2 Hd = to minimize the value of /z 2 , and require that m 2 Hu is 
positive. The minimum value of m f - c /Ms and mi /Ms for different choices of jl/Mg are gotten 
by solving m 2 c (/i G ) = and m 2 (/ic) = 0. These results are presented in Figure [7|. In this Figure 
the positivity constraints terminate at that value of which gives m 2 H = 0. 

As discussed in the above, reducing the value of m 2 Hu decreases the positivity limit on mj. . 
Consequently the fine tuning of m| with respect to m~ is also reduced. But using Equations 
|20l and |], it can be seen that decreasing m# while keeping m~,(/i G ) = and m 2 (/ic) = 
results in a larger /i, thus increasing the fine tuning with respect to //. This can also be seen 
from Figure [7[ We find, for example, that if M^p/Ms ~ 0.01, the small value fl/Ms = 0.01 
requires m^. /Ms ~ .25. For Ms = 10 TeV, this corresponds to /i ~ 100 GeV and m 4 ~. > 
2.5 TeV. A further inspection of Figure [7| shows that for the same value of M 3 ^/Ms, a value of 
mip/Ms = 0.17 is allowed (gotten by decreasing m 2 Hu ) only if jl/Ms is increased to .14. This 
corresponds to fi = 1.4 TeV for Ms = 10 TeV; this implies that A(m|; fi) ~ 930. We find that 
the limit on the initial stop masses can only be decreased at the expense of increasing fi. 

Finally, the limits become weaker if m 2 H < 0. This possibility is theoretically unattractive 
on two accounts. Firstly, a nice feature of supersymmetric extensions to the SM is that the 
dynamics of the model, through the presence of the large top quark Yukawa coupling, naturally 
leads to the breaking of the electroweak symmetry [Q. This is lost if electroweak symmetry 
14 This assumption is unnecessary. See the previous footnote. 
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Figure 5: Limits for mo/Ms from the requirement that the stop mass squareds are positive at 
the weak scale, for tan/5 = 2.2, A t $ = and assuming universal scalar masses at Mqjjt for the 
stop and Higgs scalars. The regions below the curves are excluded. 

breaking is already present at the tree-level. Secondly, the fine tuning required to obtain the 
correct Z mass is increased. From Figure [7] we infer that while reducing m 2 Hu below zero does 
reduce the limit on the initial stop masses, the value of \x increases beyond the values quoted in 
the previous paragraph, thus further increasing the fine tuning of the Z mass. This scenario is 
not discussed any further. 

We now combine the positivity analysis of this Section with the results of Section |2| to place 
lower limits on the soft scalar masses. For given values of &ll,&rr, & minimum value of Ms, 
Ms,min, is found using the results of Section ||. This is combined with the positivity analysis 
in Figure to produce the results shown in Figure |^. We also show other limits gotten by 
assuming m 2 Hu = m 2 c Q . These results are presented in Figure In Figure [lOj we also present 
the stop mass limits for different values of /i, and restrict to m% > and VSll^rr = 0.04. In 
all cases the heavy scalars were decoupled at Ms, m in, rather than 50 TeV, and so the positivity 
analysis was repeated. The value of A ti o was chosen to maximize the value of the stop masses 
at the weak scale. For completeness, the results for the cases (ri5,n w ) = (2,0) and (0,2) and 
m 2 Hu = are presented in Figure [11]. We repeat that the minimum allowable values for the stop 
masses consistent with m 2 Hu > 0, gotten by setting m 2 Hu = 0, are given in Figures || and [11]. 
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Figure 6: Limits for rrif /Ms, m^ /M s , from the requirement that the stop mass squareds are 
positive at the weak scale, for M SUS y = M GUT , tan (3 = 2.2 and assuming that rn 2 Hu = 0. The 
value of A tt0 is chosen to maximize the value of the stop soft masses at the weak scale. The 
regions below the curves are excluded. 
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Figure 7: Limits for rrif /Ms, m^c /Ms, from the requirement that the stop mass squareds are 
positive at the weak scale, for (n5,n 10 ) = (2,2), Msjjsy — Mqut, tan/5 = 2.2, and different 
values of fx/Ms- The contours end at that value of M 3t0 /Ms that gives mH u ,o/Ms = 0. The 
value of A t fl is chosen to maximize the value of the stop soft masses at the weak scale. The 
regions below the lines are excluded. 
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Figure 8: Limits for m^ and m^c , m £ - c , and from the requirement that the mass squareds 
are positive at the weak scale while suppressing Am^. It was assumed that Msusy — ^gut-, 
tan/3 = 2.2 and that m 2 H = 0. The value of At$ was chosen to maximize the value of the stop 
soft masses at the weak scale. The heavy scalar s were decoupled at the minimum value allowed 
by Amx- The regions below the lines are excluded. 

We next briefly discuss some consequences of this numerical analysis. We concentrate on the 
case n 5 = n w = 2, since this is the relevant case to consider if the supersymmetric flavour problem 
is explained by decoupling the heavy scalars. Other choices for n 5 and n w requires additional 
physics to explain the required degeneracy or alignment of any light non-third generation scalars. 
From Figures |] and || we find that for \/8l,l&rr — 0.22 and M 3j0 < 1 TeV, mf. ^7 TeV is 
required. If instead we restrict both A(m|; Mj) and A(m|; M^p) to be less than 100, then we 
must have Ms& 10 TeV and M 3i0 ^ 300 GeV. To not be excluded by Amx, we further require 
that V '5 ll<5rr^- 06. For this value of \/6llSrr = 0.06, a minimum value for mj of ~1.5— 2.5 
TeV is gotten by rescaling the results in Figures [| and || for ^/5ljJrr = 0.04 by an amount 
0.06/0.04. The range depends on the value of m H , with the lower (upper) limit corresponding 
to triHufi = ( m fc )- Thus A(m|;m~ Q ) ~ 400 — 800. This fine tuning can be reduced only 
by either increasing M^q— which increases A(m|, M^q) beyond 100— or by reducing Ms— which 
requires a smaller value for \Z5ll8rr. We conclude that unless \Z5rr5ll is naturally small, 
decoupling the heavy scalars does not provide a natural solution to the flavour problem. 
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Figure 9: Limits for mi and m t - c from the requirement that the stop mass squareds are positive 
at the weak scale while suppressing Am^. It was assumed that Msusy — Mqut, tan/3 = 2.2 
and that m 2 Hu = m~ c Q . The value of A tt0 was chosen to maximize the value of the stop soft 
masses at the weak scale. The heavy scalars were decoupled at the minimum value allowed by 
ArriK- The regions below the lines are excluded. 
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Figure 10: Limits for m ( - and m t - Ci0 from the requirement that the stop mass squareds are 
positive at the weak scale while suppressing Amjf, for (715,7110) = (2,2), \Z5ll5rr = 0.04, and 
different values of /i. The contours terminate at m 2 H Q = 0. It was assumed that M S ijsy = M GUT 
and tan/3 = 2.2. The value of A t> o was chosen to maximize the value of the stop soft masses at 
the weak scale. The heavy scalars were decoupled at the minimum value allowed by Amjf. The 
regions below the lines are excluded. 
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Figure 11: Limits for m^ , mfc from the requirement that the stop mass squareds are positive 
at the weak scale while suppressing Am^-, for the cases {115,1110) = (2,0) and (0,2). It was 
assumed that Msusy = Mqut, tan/3 = 2.2 and that m 2 Hu Q = 0. The value of A tt0 was chosen to 
maximize the value of the stop soft masses at the weak scale. The heavy scalars were decoupled 
at the minimum value allowed by Aijik- The regions below the lines are excluded. 
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To conclude this Section we discuss the constraint this analysis implies for those models which 
generate a split mass spectrum between different generations through the D-term contributions 
of the anomalous U(l) gauge symmetry [12, |16|, [OJ. In the model of set D of |15|], there are two 5s 
at 7 TeV and 6.1 TeV and two 10s at 6.1 and 4.9 TeV, respectively, so that Am K is suppressed. 
These values must be increased by a factor of 2.5 to correct for the QCD enhancement of the 
SUSY contribution to Am^, as discussed in Section |2J. To obtain a conservative bound on the 
initial stop masses from the positivity requirement, we first assume that all the heavy scalars 
have a common mass M$ = 2.5 x 5TeV= 12.5 TeV. (It would have been 5 TeV without the QCD 
correction.) Then assuming a weak scale value of the gluino mass that is less than 710 GeV 
and setting m 2 Hu 
This leads to A(mf 



Hl ), 



we find from Figure |] fl5|) that mj > 2.1 (3.6) TeV is required. 
> 580 (1700). To obtain a better bound, we repeat our analysis using 
n 5 m 2 + 3n 10 m 2 Q = ((7 TeV) 2 + (6.1 TeV) 2 + 3x(6.1 TeV) 2 + 3x(4.9 TeV) 2 ) x (2.5) 2 . It is possible 
to do this since only this combination appears in the RG analysis for (n 5 , n w ) = (2, 2). We find 



771 t,o 



(assuming rri 2 H = and the gluino mass at the weak scale is less than 710 GeV) that m ( - ~ 2.4 



TeV. In the model of [IB|, 5rr ~ 5ll ~ 0.01. To obtain a limit on the initial stop masses, we 
use the bound obtained from either Figures || or ||for Srr = 5ll ~ 0.04, and divide the limit by 
a factor of 4. By inspecting these Figures we find that this model is only weakly constrained, 
even if m 2 Hu ~ m 2 Q . We now discuss the limits in this model when 0(1) CP violating phases 
are present. To obtain the minimum value of Ms in this case, we should multiply the minimum 
value of Ms obtained from the Am^ constraint for 5ll — $rr. — 0.04 by 12.5/4; dividing by 4 
gives the result for 5ll — $rr — 0.01, and multiplying by 12.5 gives the constraint on Ms from 
e. The result is Mg^ 23 TeV. Next, we assume that M 3 is less than 300 GeV, so that the value 
of the gluino mass at the weak scale is less than 710 GeV. This gives M^^/Mg < 0.013. Using 
these values of M^p and Ms, an inspection of Figures [5] and [| implies that mi must be larger 
than 3.9 TeV to 6.9 TeV, depending on the value of r m 2 Hu . This gives A(m|;m 2 ) > 2000. In 
the model of [12], M30/M5 ~ 0.01 and rrij Q /Ms ~ 0.1. Inspecting Figures || and ^ we find 
that these values are excluded for (n 5 ,nio) = (2,2) and (0,2). The case (2,0) is marginally 
allowed. The model of [12j with (ns,nio) = (2,2) and At = was also excluded by the analysis 



of Reference [18 



5 Using Finetuning to Constrain 5 

In this section, we vary the messenger scale, Msusy, between the GUT scale and a low scale ~ 
50 TeV, and restrict the boundary values of the stop and gluino masses so that EWSB is not fine 
tuned. This gives us an upper limit to 5 if we require both positivity of the stop mass squareds at 
the weak scale and suppression of Arrix- In other words, we determine the values for (5, Msusy) 
which are allowed by the following requirements: 1. Suppression of the SUSY contribution to 
Am/f by making the mass of the first two generation scalars, M$, large. 2. Positivity of the 
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stop mass squareds and 3. Fine tuning in electroweak symmetry breaking does not exceed 1% 
or 10% (i.e., both A(m|,m~ ) and A(m 2 z , M^^) are smaller than either 100 or 10). 

An upper limit to S satisfying the above requirements is obtained as follows. For a given 
Msusy we compute, using Equations [I] and [2|, the boundary values of the stop mass, mi max , 
and the gluino mass, M^ max , such that both A(m|,m~ ) and A(m|,M 3i0 ) are equal to some 
maximum value A max which is chosen to be 100 or 10. Q Substituting these values of the bare 
stop[^] and gluino masses into the expression for the weak-scale value of the stop mass squared, 
we determine the maximum value of Ms, Ms, ma x, such that the stop mass squareds at the weak 
scale are positive. Using this value for Ms and the analysis described in Section |2.1| , an upper 
value to 5 is gotten from the Atuk constraint. This value of 5 and Msusy then satisfies the 
above-mentioned three requirements. This can be seen as follows. For the given Msusy, if 8 is 
larger than this limit, then to suppress Am^, Ms has to be larger than Ms y max- But, then to 
keep the stop mass squareds positive at the weak scale, the boundary value of either the stop or 
the gluino mass has to increase beyond rrif max or M 3tmax respectively, leading to A(m|,m~ ) or 
A(m|,M 3)0 ) larger than A max , i.e., increasing the fine tuning in EWSB. 

We show the limits on y/Sn^RR as a function of M S usy for the case (n 5 = 2, n 10 = 2) 
in Figures [12. In the top of Figure [12, m 2 Hu = is assumed. GMSB relations between the 



stop and Higgs masses are assumed in the bottom of Figure [12|. For both cases, A max = 100, 
tan/3 = 2.2 and 10 are considered. For other choices for A max , the upper limit to 5 roughly 



scales as ^A max /W0, since both m imax , M^ max and therefore M s , max scale as y/A max . 

In the case of GMSB mass relations, the boundary value of the Higgs mass and the stop 
masses are comparable for high M S usy- Since rn% u0 results in a negative contributon to the 
stop mass squared, this tends to reduce the stop mass squared at the weak scale as compared 
to the case m# = 0. Then, from the above analysis, we can see that Ms, m ax and, in turn, the 
limit on 5 is smaller for the GMSB case as compared to the case vr? H = 0. This can be seen 

15 In computing the A's, tan j3, in addition to m H (mz), should be regarded as a function of the bare parameters. 
However, this additional contribution to the A's is small for tan/3 ~ 2 and also makes the magnitude of A larger. 
We neglect this dependence which is a conservative choice. 

16 Strictly speaking, we should translate the upper bound on to? Q into an upper bound on fh~ using to? Q = 
m? q+cdY^.Z^.o+^Cdj *- e -i t° that combination appearing in the positivity constraint. Instead, we use the same 
bound for both to? Q and m? Q . This is reasonable, since cd is generally small (^0.05), and Dy.o ~ 0(m 2 ). In any 
case, this effect is in the opposite direction for t and t c . In the case that £ D ^ 0, a slightly larger (O(30%)) value 
for 6 may be allowed as compared to £o = 0. This is because if Cd < 0, the maximum value for m~ c Q is larger than 
w ? max' This, in turn, allows for a larger value of Ms, and hence 8. Naturalness considerations limit \Cd\, though. 
The EWSB relation for m 2 z , Equation ^, contains a term linear in C,d- Requiring that A(w|, Cd) < 100 implies 
that |£d |^C-D, max = (900 GeV) 2 . Thus for a high scale of supersymmetry breaking, the upper bound on m? c 
may be increased to m| Q ~ ™\ max + j(D,max ~ \ m \ max -> while maintaining A(m|, m? c f) ) = A(m|, Cd) = 100. 
This roughly translates into an increase of ~ ^/5/3 = 1-3 in the limit to S. The actual limit will be smaller, 
since with this choice of sign for Cd, the positivity constraint for the left-handed stop is now stronger. It is thus 
reasonable to require that the maximum value of to? _ be comparable to to? 

1 r^,0 ZA^max 
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by comparing the top and bottom of Figure |12[ 

In Figure [13| the limits on Srr and 5ll for (n 5 = 2, n w = 0) and (n 5 = 0, nio = 2) are shown, 
respectively. We assume m 2 H = and consider tan/? = 2.2 and 10. If we choose A max to 
be 100, then we get a constraint on 5 (5 ~ 0.5) only for high values of M S usy- So, we choose 
instead A max to be 10. 

We have checked that, for tan/3 = 10, the limits on the boundary value of the stop mass 
from requring positivity of the mass squared at the weak scale do not differ by more than a few 
percent from the case tan/5 = 2.2 (for the same values of the gluino and heavy scalar masses). 
However, the fine tuning of EWSB for the same gluino and stop mass is smaller for tan/3 = 10 
as compared to tan/3 = 2.2. This is because, for tan/3 = 10, Xt is smaller than in the case 
tan/3 = 2.2. Hence the sensitivity of the weak scale value of rn 2 Hu to m~ and M 3; o is smaller. 
Also, the tan 2 /3/(tan 2 f3 — 1) factor in Equation |2| is smaller, further reducing the sensitivity 
of m 2 z to m 2 Q and M 3j0 . In other words, for tan/3 = 10, m^ max and M^ max are larger so that 
Ms max and, in turn, the limit on 8 is larger. This can be seen in Figures O and [13. 



6 Conclusions 

In this paper we have studied whether the SUSY flavor problem can be solved by making 
the scalars of the first and second generations heavy, with masses Ms (~few TeV), without 
destabilising the weak scale. If the scale, Msusy, at which SUSY breaking is mediated to the 
SM scalars is close to the GUT scale, then the heavy scalars drive the light scalar (in particular 
the stop) mass squareds negative through two-loop RG evolution. In order to keep the mass 
squareds at the weak scale positive, the initial value of the stop (and other light scalar) soft 
masses, mj. , must typically be ~ 1 TeV, leading to fine tuning in EWSB. We included two new 
effects in this analysis: the effect of \ t in the RGEs which makes the stop mass squareds at the 
weak scale more negative and hence makes the constraint on the initial value stronger, and the 
QCD corrections to the SUSY box diagrams which contribute to K — K mixing. 

Some results of our analysis for Msusy = Mqut can be summarized as follows. We restrict 
the gluino mass (at the weak scale) to be less than about 710 GeV, so that the fine tuning of m 2 z 
with respect to the bare gluino mass, M 3j0 , is not worse than 1%. This requires that M 30 ^300 
GeV. We also assume that m 2 Hu = to maximize the value of the stop masses at the weak 
scale. We find that if V^ll^rj? = 0.22 then M s > 40 TeV is required to be consistent with 
ArriK- With these assumptions, this implies that for M 3j0 less than 1 TeV, m 4 ~ > 6.5 TeV is 
needed to not break colour and charge at the weak scale. Even for \^5ll^rr = 0.04, we find 
that we need Ms ~ 7 TeV. This implies that m t - > 1 TeV is required if M 3 < 300 GeV. This 
results in a fine tuning of ~ 1%. For 5ll = 1 an d 8rr = 0, we find that M s ~ 30 TeV and 
m ifl > 4.5 TeV. For 5 LL = 0.22 and 5 RR = 0, we find that M s ~ 7 TeV and m i0 > 1 TeV. For 
5ll = and 5rr = 1, we find that Ms ~ 30 TeV and mf c0 > 2.5 TeV. The constraints are 
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Figure 12: Maximum value for (Sll^rr) 1 ^ 2 that is consistent with A(m|, M^p) < 100, 
A(m|,m~ ) < 100 and (ArriK)susY < (Amif)^. Two boundary conditions are considered: 
m 2 H o = (top) and gauge- mediated relations (bottom). Two values for tan/5 are considered. 
The value of A tj o was chosen to maximize the value of the stop masses at the weak scale. 



32 



to 



to 




d 1 



10 



A(M )=10 + A(m,)=10 



10 



15 



20 



25 



10 



15 



20 



25 



n 5 =2, n w =0 

• tan $=2.2 
T tan $=10. 




30 



35 



In M„,,„ v /m t 

SUSY top 





1 1 1 1 1 1 1 1 1 1 1 1 1 

A(M )=10 + A(m t )=10 


1 1 1 1 1 1 

n 5 =0, n w =2 I 






• tan $=2.2 






T tan $=10. 


1 1 


i i i i i i i i i i i i i i 





30 



35 



ltlM SUSY /m top 



2 ~ ) 

tx>) 



< 



Figure 13: Maximum value for 5ll, $rr that is consistent with A(m|, M 3 ) < 10, A(m|, m 
10 and (ArriK)susY < (Am^ezp- It was assumed that m 2 Hu0 = 0. Two values for tan/? are 
considered. The value of A t>0 was chosen to maximize the value of the stop masses at the weak 
scale. 
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weaker for smaller values of 5. In a realistic model, m 2 Hu might be comparable to m~ and the 
constraints on mj in this case are stronger. This is also discussed. We note that independent 
of the constraint from K — K mixing, our analysis can be used to check the phenomenological 
viability of any model that has heavy scalars. We also discuss the phenomenological viability of 
the anomalous D— term solution, and find it to be problematic. 

We then considered the possibility that M SUSY = M s . In this case, there is no RG log 
enhancement of the negative contribution of the heavy scalar masses to the light scalar masses. 
For this case, we computed the finite parts of the two-loop diagrams and used these results as 
estimates of the two-loop contribution of the heavy scalars to the light scalar soft mass squareds. 
We then combined these results with the constraints from K — K mixing to obtain lower limits on 
the boundary values of the stops. As an example, we assumed gauge mediated SUSY breaking 
boundary conditions for the light scalars. If n 5 7^ riio then one of the selectron masses, rather 
than the stop masses, provides the stronger constraint on m t ~. once gauge-mediated boundary 
conditions are used to relate mgc j0 and mi to m~ ti . Some of our results can be summarized as 
follows. We restrict the gluino mass at the weak scale to be less than about 2.3 TeV, again to 
avoid more than 1% fine tuning of m\ with respect to the gluino mass. For \/5Zl7rr — -22 we 
find that mj. > 1.4 TeV is required. The fine tuning of m 2 z with respect to the stop mass is 
~ 1.5% in this case. For the cases 5ll = and Srr = 1, and 5ll = 1 and 5rr = we find that 
mi o ~ 1 TeV. As before, the constraints on mj for smaller values of 5 are weaker than ~ 1 TeV. 
Again, we emphasize that the constraints in an actual model of this low-energy supersymmetry 
breaking scenario could be different, and our results should be treated as estimates only. We 
also briefly discuss the CP violating constraints from e, and find that these limits increase by a 
factor of ~ 12 if 0(1) phases are present. 

Finally, in Section 5 the scale of supersymmetry breaking is varied between 50 TeV and 
2 x 10 16 GeV. Uppers bounds to 5, that are consistent with positivity of the light scalar masses, 
naturalness in electroweak symmetry breaking, and (Am^)^, are obtained. These results are 
summarized in Figures [12| and O. 
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8 Appendix: Two-loop calculation 

In this Appendix we discuss the two-loop contribution of the heavy scalar soft masses to the 
light scalar soft masses. These contributions can be divided into two classes. In the first class, 
a vev for the hypercharge D-term is generated at two-loops. The Feynman diagrams for these 
contributions are given in Figure and are clearly ~ a.\cti. These diagrams are computed in 
a later portion of this Appendix. In the other class, the two- loop diagrams are ~ af. These 



have been computed by Poppitz and Trivedi [30]. So, we will not give details of this computation 



which can be found in their paper. However, our result for the finite parts of these diagrams 
differs slightly from theirs and we discuss the reason for the discrepancy. When one regulates 



the theory using dimensional reduction [£5], |26fl (compactifying to D < 4 dimensions), the vector 
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field decomposes into a D-dimensional vector and 4 — D scalars, called e-scalars, in the adjoint 
representation of the gauge group. Thus the number of Bose and Fermi degrees of freedom in 
the vector multiplet remain equal. The e-scalars receive, at one-loop, a divergent contribution 
to their mass, proportional to the supertrace of the mass matrix of the matter fields. Neglecting 
the fermion masses, this contribution is 

a 2 

5m 2 e = -—(- + ln4vr - 7 )(n 5 + 3n w )M 2 s . (25) 
47T e 

In our notation D = 4 — e. Poppitz and Trivedi choose the counterterm to cancel this divergence 
in the MS scheme, i.e., the counterterm consists only of the divergent part, proportional to 1/e. 
When this counterterm is inserted in a one-loop e-scalar graph with SM fields (scalars) as the 
external lines , one obtains a divergent contribution to the SM scalar soft masses (the 1/e of the 
counterterm is cancelled after summing over the e adjoint scalars running in the loop). Poppitz 
and Trivedi use a cut-off, A[/y, to regulate this graph, giving a contribution from this graph that 
is: 

m 2 = -5>6 + 3n 10 )C^(^Ml\nA 2 uv (26) 

with no finite part. We, on the other hand, choose the e-scalar mass counterterm in the MS 
scheme, i.e., proportional to 2/e — 7 + ln47r (where 7 ~ 0.58 is the Euler constant) and use 
dimensional reduction to regulate the graph with the insertion of the counterterm. This gives a 
contribution 

< = -EK + 3n 10 )C^(^) 2 M|(--7 + ln47r) 2 e 
^ 16 V 7r / e 

= -EK + 3n 10 )C^^) 2 M|(2/e-27 + 21n47r) (27) 

In the first line the first factor of (2/e — 7 + ln47r) is from the counter-term insertion, the second 
factor is the result of the loop integral, and the over-all factor of e counts the number of e-scalars 
running in the loop. In the MS scheme, i.e., after subtracting 2/e — 7 + ln47r, we are left with 
a finite part[] proportional to — 7 + ln47r. The remaining diagrams together give a finite result 
and we agree with Poppitz and Trivedi on this computation. Our result for the finite part of the 
two-loop diagrams (neglecting the fermion masses) is 

ml fini M = 4( ln(4ff) ~ 7 + y ~ 2 ~^7?)) 

xE(^)V 5 + 3n 10 )QM| (28) 

whereas the Poppitz- Trivedi result does not have the ln(47r) — 7 in the above result. The com- 
putation of the two-loop hypercharge D-term, which gives contribution to the soft scalar mass 
17 The same finite part is obtained in the MS scheme, regulated with DR' . 
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Figure 14: Mixed two-loop corrections to the scalar mass. Wavy lines, wavy lines with a straight 
line through them, solid lines, and dashed lines denote gauge boson, gaugino, fermion and scalar 
propagators,respectively. The double-line denotes the hypercharge D-term propagator. 

squareds proportional to a\a s and a\a 2 (i.e., the "mixed" two-loop contributon) is discussed 
below in detail. 

Two-loop hypercharge D-term 



We compute the two-loop diagrams of Figure [L4| in the Feynman gauge and set all fermion 
and gaugino masses to zero. It is convienent to calculate in this gauge because both the scalar 
self-energy and the Dy-term vertex corrections are finite at one-loop and thus require no counter- 
terms. We have also computed the two-loop diagrams in the Landau gauge and have found that 
it agrees with the calculation in the Feynman gauge. The calculation in the Landau gauge 
requires counter-terms, is more involved, and hence the discussion is not included. Finally, 
in the calculation a global £77(5) symmetry is assumed so that a hypercharge JJ-term is not 
generated at one-loop |L7], [13]]. 



The sum of the four Feynman diagrams in Figure [L4] is given in the Feynman gauge by 



-if[ D>f = dg'YfTrY^glC^him*) - 4J 2 (m?) + h{m 2 )) (29) 

A,i 

where the trace is over the gauge and flavour states of the particles in the loops. If the particles 
in the loop form complete 5 and 10 representations with a common mass Ms, the sum simplifies 
to 

- iti DJ = i^axY f (n 5 - n 10 )(^a 3 - ^a 2 - -^ ai )(4h(M 2 s ) - 4J 2 (M|) + h(M 2 s )). (30) 
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The functions h, I 2 and J 3 are 

2 r d D p r d D k 1 (2p-k) 2 1 

l[m ] J {2n) D J {2tt) d (p 2 - m 2 ) 2 k 2 (p - k) 2 - m 2 ' 1 ' 

T( 2^ f d D p r d D k 1 k 2 -k-p 1 

2lm j "~ J (2n) D J (2n) D (p 2 - m 2 ) 2 k 2 (p - k) 2 ' { ' 

T( 2X /• 1 f dD( l 1 /OOA 

3lm } " 7 (2tt) d {k 2 - m 2 ) 2 J (2n) D q 2 - m 2 ' { ' 
We now compute these functions. 
Evaluating h 

After a Feynman parameterization and performing a change of variables, h = J\ + h, where 

Jl(m2) = r(3) f ^ - X) I^J ^Y 2 {p 2 -{m 2 -x{l-x)k 2 )Y (M) 

and 

*^ - r(3) / V - ,)(2, - If / ^ / ^- D {jf _ (m2 _l (i _ x)k2)r (35) 
After some algebra we find that 

Mm*) = I ^^(rnT- 3 ^- [ B(2 - D/2, 3 - D/2), (36) 

J 2 (m 2 ) = T ^~^\ m 2 ) D ^x{AB^-D/2,2-D/2)-AB{2-D/2,2-D/ 

(37) 

where q) = Y[p]T[q]/Y[p + g] is the usual Beta function. 
Combining these two results gives 



1^) = £ |^K)^ 3 ^5(3 " D/2, 2 - D/2). (38) 



Evaluating h 



d D P r d D k 1 k 2 -k- P 1 



Evaluating 7 3 



(47T) 



(27t)^7 (2tt) d (p2 - m 2)2 p ( p _A;)2 
' r(3- J D)(m 2 ) z? - 3 B(D/2,l- J D/2). 



y tFfc 1 y d D q 1 
7 (2n) D (k 2 - m 2 ) 2 J (2n) D q 2 -m 2 
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We may now combine I i: I 2 and I 3 to obtain 

T(m 2 ) = 4/i (m 2 ) -4/ 2 (m 2 ) + him 2 ) 
(m 2 ) D ~ 3 ,/1-D 



~ X{A {d^2 B{3 ~ D/2 > 2 ~ D/2) ~ B{D/2 > 1 " D/2) ) r(3 ~ D) 



(An) 

1 -T(2-D/2) 2 ). 



D/2-1 

Writing D = 4 — e and expanding in e gives 

T(m2) = (16^p (l + ( 6 " ^ 2 + 4 ( ln ^ 47r ) - 7) - 4 him 2 ) m 2 + 0(e)) . (39) 



In the MS scheme the combination 2/e + ln(47r) — 7 is subtracted out. The finite piece that 
remains is 



(16tt 2 ) 



f6 - -7t 2 + 2(ln(47r) - 7) - 41nm 2 ) m 2 . (40) 

7T 2 Y V 3 / 



Thus in the MS scheme 



5(16tt 2 ) 2 

.3 1 / 2 o _ ,M; 



2 



= «5^«i(^)(n 5 - n 10 )y/ (6 - ^tt 2 + 2(1h(4tt) - 7) - 41n(-^)J 

x (gO&W - ^2^) - ^aiWj M 2 (41) 
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